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O ■ Abstract 
O; 

The multiple Doppler readouts available on the Laser Interferometer Space Antenna (LISA) 
permit simultaneous formation of several interferometric observables. All these observables are 

Hi ! 

■^j- ■ independent of laser frequency fluctuations and have different couplings to gravitational waves and 

^NJ . to the various LISA instrumental noises. Within the functional space of interferometric combina- 

> ; 

>> \ ' tions LISA will be able to synthesize, we have identified a triplet of interferometric combinations 

m : 

© ; that show optimally combined sensitivity. As an application of the method, we computed the sen- 

sitivity improvement for sinusoidal sources in the nominal, equal-arm LISA configuration. In the 
part of the Fourier band where the period of the wave is longer than the typical light travel-time 

O 

across LISA, the sensitivity gain over a single Michelson interferometer is equal to \/2. In the 
mid-band region, where the LISA Michelson combination has its best sensitivity, the improvement 



over the Michelson sensitivity is slightly better than y/2, and the frequency band of best sensitivity 



is broadened. For frequencies greater than the reciprocal of the light travel-time, the sensitivity 
improvement is oscillatory and ~ v3, but can be greater than \/3 near frequencies that are integer 
multiples of the inverse of the one-way light travel-time in the LISA arm. 
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I. INTRODUCTION 



The Laser Interferometer Space Antenna (LISA) is a deep-space mission, jointly proposed 
to the National Aeronautics and Space Administration (NASA) and the European Space 
Agency (ESA), to detect and study gravitational radiation in the millihertz frequency band. 

LISA will use coherent laser beams exchanged between three remote, widely separated, 
spacecraft. Modeling each spacecraft as carrying lasers, beam splitters, photo-detectors, and 
drag-free proof masses on each of two optical benches, it has been shown jlj], Q], that the 
six measured time series of Doppler shifts of the one-way laser beams between spacecraft 
pairs, and the six measured shifts between adjacent optical benches on each spacecraft, can 
be combined, with suitable time delays, to cancel the otherwise overwhelming frequency fluc- 
tuations of the lasers (Au/u ~ 10 -13 /\/Hz), and the noise due to the mechanical vibrations 
of the optical benches (which could be as large as Au/u ~ 10 _16 /v / Hz). The achievable 
strain sensitivity level h ~ 10~ 21 / a/Hz is set by the buffeting of the drag-free proof masses 
inside each optical bench, and by the shot noise at the photo-detectors. 

In contrast to Earth-based, equi-arm interferometers for gravitational wave detection, 
LISA will have multiple readouts. The data they generate, when properly time shifted 
and linearly combined, provide observables that are not only insensitive to laser frequency 
fluctuations and optical bench motions, but also show different couplings to gravitational 
radiation and to the remaining system noises , 0| . The space of all possible inter- 

ferometric combinations can be generated by properly combining four generators and 
it has been shown to be an algebraic module [5j. In this paper we derive from first princi- 
ples a particular combination of these generators which, for a given waveform and source 
location in the sky, give maximum signal-to-noise ratio. In this respect, LISA should no 
longer be regarded as a single detector, but rather as an array of interferometers working in 
coincidence. 

An outline of the paper is presented here. Section II provides a brief summary of Time- 
Delay Interferometry, the data processing technique needed to remove the frequency fluctu- 
ations of the six lasers used by LISA. After showing that the entire set of interferometric 
combinations can be derived by properly combining four generators, (a,/?, 7, £), in Section 
III we turn to the problem of optimization of the Signal-to-Noise ratio (SNR) within this 
functional space. As an application, we apply our results to the case of sinusoidal signals 
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randomly polarized and randomly distributed on the celestial sphere. We find that the 
standard LISA sensitivity figure derived for a single Michelson Interferometer 2] can be 
improved by a factor of v2 in the low-part of the frequency band, and by more than y/3 
in the remaining part of the accessible band. In Section IV we present our comments and 
conclusions. 



II. TIME-DELAY INTERFEROMETRY FOR LISA 

Figure 1 shows the overall LISA geometry. The spacecraft are labeled 1, 2, 3 and distances 
between pairs of spacecraft are L\, L 2 , L 3 , with Li being opposite spacecraft i. Unit vectors 
between spacecraft are "hi, oriented as indicated in Figure 1. We similarly index the relative 
frequency fluctuations j/y data to be analyzed: y^% is the relative frequency fluctuations 
time series measured at reception at spacecraft 1 with transmission from spacecraft 2 (along 
L 3 ). Similarly, y 2 % is the corresponding time series derived from reception at spacecraft 1 
with transmission from spacecraft 3. The other four one-way relative frequency time series 
from signals exchanged between the spacecraft are obtained by cyclic permutation of the 
indices: 1 — ► 2 — > 3 — ► 1. The useful notation for delayed data streams will be also used: 
1/31,2 = V3i{t - L 2 ), ?/3i,23 = V3i{t - L 2 - L 3 ) = #31,33, etc. (units in which c = 1). Six 
more Doppler difference series result from laser beams exchanged between adjacent optical 
benches within each spacecraft; these are similarly indexed as z^ = 1,2,3 ; i ^ j) (see 
^ and br details). 

The light paths for the y^s and z^s can be traced in Figure 2. An outgoing light 
beam transmitted to a distant spacecraft is routed from the laser on the local optical bench 
using mirrors and beam splitters; this beam does not interact with the local proof mass. 
Conversely, an incoming light beam from a distant spacecraft is bounced off the local proof 
mass before being reflected onto the photo-detector where it is mixed with light from the 
laser on that same optical bench. These data are denoted y%i and y 2 i in Figure 2. Beams 
exchanged between adjacent optical benches however do precisely the opposite. Light to be 
transmitted from the laser on an optical bench is first bounced off the proof mass it encloses 
and then directed via fiber optics to the other optical bench. Upon reception it does not 
interact with the proof mass there, but is directly mixed with local laser light. They are z 3 i 
and z 2 \ in Figure 2. 
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FIG. 1: Schematic LISA configuration. Each spacecraft is equidistant from the point O, in the plane 
of the spacecraft. Unit vectors hi point between spacecraft pairs with the indicated orientation. 
At each vertex spacecraft there are two optical benches (denoted 1, 1*, etc.), as indicated. 



er observables yij, Zij with specific time signatures. 
3) Hi) an d we refer the reader to those papers for a 



The frequency fluctuations introduced by the gravitational wave signal, the lasers, the 
optical benches, the proof masses, the fiber optics, and the measurement itself at the photo- 
detector (shot noise) enter into the Dopp 
They have been derived in the literature 
detailed discussion. The Doppler data y^, are the fundamental measurements needed to 
synthesize all the interferometric observables unaffected by laser and optical bench noises. 

The simplest such combination, the totally symmetrized Sagnac response (, uses all the 
data of Figure 2 symmetrically 

C = 1/32,2 - 1/23,3 + 1/13,3 ~ 1/31,1 + 1/21,1 ~ 1/12,2 

I, . 
+ — v — z 13,21 + ^23,12 ~ 221,23 + £31,23 — ^32, 13 + Z 12,13J 

~\~2^~ Z32 ' 2 ~^ Zl2,2 ~ Zl3 ' 3 Z23 ' 3 ~~ 221,1 Zs1 '^ ' W 

and its transfer functions to gravitational waves and instrumental noises were derived in Q], 
and |2| respectively. In particular, ( has a "six-pulse response" to gravitational radiation, 
i.e. a 5-function gravitational wave signal produces six distinct pulses in ( which are 
located with relative times depending on the arrival direction of the wave and the detector 
configuration. 

Together with (, three more interferometric combinations, (a,(3,j), jointly generate the 
entire space of interferometric combinations [lj], Their expressions in terms of the 
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FIG. 2: Schematic diagram of proof-masses-plus-optical-benches for a LISA spacecraft. The left- 
hand bench reads out the Doppler signals y^i and Z31. The right hand bench analogously reads out 
■1/21 and Z2i- The random velocities of the two proof masses and two optical benches are indicated 
(lower case Vi for the proof masses, upper case Vi for the optical benches.) 

measurements yij, Zij are as follows 

« = 2/21 — 2/31 + 2/13,2 - 2/12,3 + 2/32,12 ~ 2/23,13 

1, . 
— 2 v Zl3 > 2 + ^13,13 + Z21 + ^21,123 + ^32,3 + ^32,12 J 

+ 2"( 2 23,2 + ^23,13 + ^31 + ^31,123 + ^12,3 + £12,12) , (2) 

with (3, and 7 derived by permuting the spacecraft indices in a. Like in the case of (, a 
5-function gravitational wave produces six pulses in a, j3, and 7. 

We should remind the reader that the four interferometric responses (a, /3, 7, () satisfy 
the following relationship 

C - C123 = a 1 - a 23 + (3,2 - (3,i3 + 7,3 - 7,12 • (3) 

Jointly they also give the expressions of the interferometric combinations derived in 

the Unequal-arm Michelson (X, Y, Z), the Beacon (P, Q, R), the Monitor (E, F, G), and the 
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Relay (U, V, W) responses 



Xi 


= a 23 - (3,2 - 7,3 + C , 


(4) 


p 


= C - «, i , 


(5) 


E 


= ol - C,i , 


(6) 


U 


= 1,1-13, 


(7) 



with the remaining expressions obtained from equations (JU El EH EI) by permutation of the 
spacecraft indices. All these interferometric combinations have been shown to add robustness 
to the mission with respect to failures of subsystems, and potential design, implementation, 
or cost advantages 

HQ. 



III. OPTIMAL SENSITIVITY FOR LISA 



All the above interferometric combinations have been shown to individually have rather 
different sensitivities |2(|, as a consequence of their different responses to gravitational radia- 
tion and system noises. LISA has the capability of simultaneously observing a gravitational 
wave signal with many different interferometric combinations, all having different antenna 
patterns and noises. We should thus no longer regard LISA as a single detector system but 
rather as an array of gravitational wave detectors working in coincidence. This suggests 
that the presently adopted LISA sensitivity could be improved by optimally combining the 
four generators (a,(3,j, £). In mathematical terms this can be restated by saying that we 
should be able to find that particular combination of the four generators that has maximum 
signal-to-noise ratio to a given gravitational wave signal. 

In order to proceed with this idea, let us consider the following linear combination rj(f) 
of the Fourier transforms of the four generators (a, (3, 7, () 

v (f) = ai (f,X) a(f) + a 2 (/,A) Mf) + as(/,X) i(f) + a 4 (f,X) ((f) , (8) 

where the {a«(/, A)}^ =1 are arbitrary complex functions of the Fourier frequency /, and of a 
vector A containing parameters characterizing the gravitational wave signal (source location 
in the sky, waveform parameters, etc.) and the noises affecting the four responses (noise 
levels, their correlations, etc.). For a given choice of the four functions {cij}f =1 , i] gives an 
element of the functional space of interferometric combinations generated by (a, (3, 7, C)- 
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Our goal is therefore to identify, for a given gravitational wave signal, the four functions 
{ai}j =l that maximize the signal-to-noise ratio, SNR 2 , of the combination rj 



fu 



\a x a s + a 2 P s + a 3 7 S + a 4 ( s 



SNR 2 = / ^ ' z 1 s ' ^ i s ' ^ df . (9) 
■Afi «„ + a 2 /3 n + a 3 7„ + a 4 Cn| 2 ) 

In equation (JUJ) the subscripts s and n refer to the signal and the noise parts of (5, (3, 7, £) 
respectively, the angle brackets represent noise ensemble averages, and the interval of inte- 
gration (fi, f u ) corresponds to the frequency band accessible by LISA. 

Before proceeding with the maximization of the SNR 2 we may notice from equation (J3J) 
that the Fourier transform of the totally symmetric Sagnac combination, (, multiplied by 
the transfer function l — e 2m ^ Ll+L2+L ^ can be written as a linear combination of the Fourier 
transforms of the remaining three generators (a, (3, 7). Since the signal-to-noise ratio of 77 
and (1 — e 2m ^ Ll+L2+L '^)r] are equal, we may conclude that the optimization of the signal- 
to-noise ratio of r] can be performed only on the three observables a, /3, 7. This implies the 
following redefined expression for SNR 2 

SNR 2 V = t |ai 5 * + a 4 + a 3 7,| 2 df . (10 ) 

The SNR 2 can be regarded as a functional over the space of the three complex functions 
{ a i}i=i, an d the particular set of complex functions that extremize it can of course be derived 
by solving the associated set of Euler-Lagrange equations. 

In order to make the derivation of the optimal SNR easier, let us first denote by x^ s ^ 
and x^™) the two vectors of the signals (5 S , /3 S , 7,) and the noises (a„, /3 n >7n) respectively. 
Let us also define a to be the vector of the three functions {oj}f =1 , and denote with C the 
hermitian, non-singular, correlation matrix of the vector random process x n 

(C) rt = (x^x^*) . (11) 

If we finally define (A)y to be the components of the hermitian matrix x-^x^*, we can 
rewrite SNR 2 in the following form 

snr '= rsirS*' < i2) 

J fl A r^rt A t 

where we have adopted the usual convention of summation over repeated indices. Since the 
noise correlation matrix C is non-singular, and the integrand is positive definite or null, 



the stationary values of the signal-to-noise ratio will be attained at the stationary values of 
the integrand, which are given by solving the following set of equations (and their complex 
conjugated expressions) 





, k =1,2, 3. (13) 



After taking the partial derivatives, equation (|13|) can be rewritten in the following form 



3p Apg3g 



(a*), , * = 1,2,3 (14) 



which tells us that the stationary values of the signal-to-noise ratio of r\ are equal to the 
eigenvalues of the the matrix C _1 ■ A. The result in equation ([T3 l is well known in the 
theory of quadratic forms, and it is called the Rayleigh's principle^ ] , 0. 

In order now to identify the eigenvalues of the matrix C _1 • A, we first notice that the 
3x3 matrix A has rank 1. This implies that the matrix C -1 ■ A has also rank 1, as it is 
easy to verify. Therefore two of its three eigenvalues are equal to zero, while the remaining 
non-zero eigenvalue represents the solution we are looking for. 

The analytic expression of the third eigenvalue can be obtained by using the property 
that the trace of the 3x3 matrix C _1 • A is equal to the sum of its three eigenvalues, and 
in our case to the eigenvalue we are looking for. From these considerations we derive the 
following expression for the optimized signal-to-noise ratio SNR^ o ^ t 

SNRl pt = (C-%- xW df. (15) 



We can summarize the results derived in this section, which are given by equations (|10|15j) . 
in the following way: 

(i) among all possible interferometric combinations LISA will be able to synthesize with its 
four generators a, /3,7, £, the particular combination giving maximum signal-to-noise ratio 
can be obtained by using only three of them, namely (a, (3, 7); 

(ii) the expression of the optimal signal-to-noise ratio given by equation (|TH|) implies that 
LISA should be regarded as a network of three interferometer detectors of gravitational 
radiation (of responses (a, (3, 7)) working in coincidence Q]. 



S 



A. Application 



As an application of equation (|15p. here we calculate the sensitivity that LISA can reach 
when observing sinusoidal signals uniformly distributed on the celestial sphere and of random 
polarization. In order to calculate the optimal signal-to-noise ratio we will also need to use a 
specific expression for the noise correlation matrix C. As a simplification, we will assume the 
LISA arm-lengths to be equal to its nominal value L = 16.67 sec, the optical-path noises to 
be equal and uncorrelated to each other, and finally the noises due to the proof-mass noises 
to be also equal, uncorrelated to each other and to the optical-path noises. Under these 
assumptions the correlation matrix becomes real, its three diagonal elements are equal, and 
all the off-diagonal terms are equal to each other, as it is easy to verify by direct calculation 
P). The noise correlation matrix C is therefore uniquely identified by two real functions, S a 
and S a /3, in the following way 

y Sa/3 S a /3 S a J 

The expression of the optimal signal-to-noise ratio assumes a rather simple form if we 
diagonalize this correlation matrix by properly "choosing a new basis" . There exists an or- 
thogonal transformation of the generators (a, (3, 7) which will transform the optimal signal- 
to-noise ratio into the sum of the signal-to-noise ratios of the "transformed" three inter- 
ferometric combinations. The expressions of the three eigenvalues {//j}? =1 (which are real) 
of the noise correlation matrix C can easily be found by using the algebraic manipulator 
Mathematica jl(|, and they are equal to 

^1 = ^2 = S a — S af 3 , /i 3 = S a + 2 S a p . (16) 

Note that two of the three real eigenvalues, (/ii,yU 2 ), are equal. This implies that the eigen- 
vector associated to /X3 is orthogonal to the two-dimensional space generated by the eigen- 
value fix, while any chosen pair of eigenvectors corresponding to fii will not necessarily be 
orthogonal. This inconvenience can be avoided by choosing an arbitrary set of vectors in 
this two-dimensional space, and by ortho-normalizing them. After some simple algebra, we 
have derived the following three ortho-normalized eigenvectors 

vi = -^ (-1,0,1) , v 2 = i= (1,-2,1) , v 3 = -^= (1,1,1). (17) 



Equation (J 17)) implies the following three linear combinations of the generators (a, (3, 7) 

A=-^(7-a) , E = ^ (5-2^+7) , T=i= (5 + ^ + 7), (18) 

where A, E, and T are italicized to indicate that these are "orthogonal modes" . Although the 
expressions for the modes A and E depend on our particular choice for the two eigenvectors 
( v i, v 2), it is clear from our earlier considerations that the value of the optimal signal-to- 
noise ratio is unaffected by such a choice. From equation ()18|) it is also easy to verify that 
the noise correlation matrix of these three combinations is diagonal, and that its non-zero 
elements are indeed equal to the eigenvalues given in equation ()16J) . 

In order to calculate the sensitivity corresponding to the expression of the optimal signal- 
to-noise ratio, we have proceeded similarly to what was done in , and described 
in more detail in We assume an equal-arm LISA (L = 16.67 light seconds), and 
take the one-sided spectra of proof mass and aggregate optical-path-noises (on a single 



3, 0; 
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link), expressed as fractional frequency fluctuation spectra, to be 
2.5 x 1(T 48 [f/lHzY 2 Hz' 1 and S° pUcal path = 1.8 x 1(T 37 [f/lHzf Hz~\ respectively. We 
also assume that aggregate optical path noise has the same transfer function as shot noise. 

The optimum SNR is the square root of the sum of the squares of the SNRs of the three 
"orthogonal modes" (A, E, T). To compare with previous sensitivity curves of a single LISA 
Michelson interferometer, we construct the SNRs as a function of Fourier frequency for 
sinusoidal waves from sources uniformly distributed on the celestial sphere. To produce the 
SNR of each of the (A, E, T) modes we need the gravitational wave response and the noise 
response as a function of Fourier frequency. We build up the gravitational wave responses 
of the three modes (A,E,T) from the gravitational wave responses of (a, /3, 7). For 7000 
Fourier frequencies in the ~10 -4 Hz to ~1 Hz LISA band, we produce the Fourier transforms 
of the gravitational wave response of (a, [3, 7) from the formulas in [lj], [llj]. The averaging 
over source directions (uniformly distributed on the celestial sphere) and polarization states 
(uniformly distributed on the Poincare sphere) is performed via a Monte Carlo method. 
From the Fourier transforms of the (a, /3, 7) responses at each frequency, we construct the 
Fourier transforms of (A, E, T). We then square and average to compute the mean-squared 
responses of (A, E, T) at that frequency from 10 4 realizations of (source position, polarization 
state) pairs. 

The noise spectra of (A, E, T) are determined from the raw spectra of proof-mass and 
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optical-path noises^and the transfer functions of these noises to (A, E, T). Using the transfer 
functions given in |2j, the resulting spectra are equal to 

SaU) = SeU) = 16 sin 2 (7r/L) [3 + 2cos(27r/L) + cos(47r/L)]Sf oof mass (f) 

+ 8 sin 2 (7r/L) [2 + cos(2tt/L)] S^ tical path (f) , (19) 



S T (f) = 2[l + 2cos(2vr/L)] 2 [4 sm 2 (n fL)S^°° f mass + S° ptical path (f)} . (20) 

Let the amplitude of the sinusoidal gravitational wave be h. The SNR for, e.g. A, SNRa, 
at each frequency / is equal to h times the ratio of the root-mean-squared gravitational wave 
response at that frequency divided by yS~A~(f) B, where B is the bandwidth conventionally 
taken to be equal to 1 cycle per year. Finally, if we take the reciprocal of SNRa/H and 
multiply it by 5 to get the conventional SNR = 5 sensitivity criterion, we obtain the 
sensitivity curve for this combination which can then be compared against the corresponding 
sensitivity curve for the equal-arm Michelson Interferometer. 

In Figure 3 we show the sensitivity curve for the LISA equal-arm Michelson response 
(SNR = 5) as a function of the Fourier frequency, and the sensitivity curve from the 
optimum weig hting of the data described above: hh/y/SNR\ + SNR% + SNR 2 T . The SNRs 
were computed for a bandwidth of 1 cycle/year. Note that at frequencies where the LISA 
Michelson combination has best sensitivity, the improvement in signal-to-noise ratio provided 
by the optimal combination is slightly larger than y2. 

In Figure 4 we plot the ratio between the optimal SNR and the SNR of a single Michelson 
interferometer. In the long- wavelength limit, the SNR improvement is y2. For Fourier fre- 
quencies greater than or about equal to l/L, the SNR improvement is larger and varies with 
the frequency, showing an average value of about \/3- In particular, for bands of frequencies 
centered on integer multiples of l/L, SNR T contributes strongly and the aggregate SNR in 
these bands can be greater than 2. In order to better understand the contribution from the 
three different combinations to the optimal combination of the three generators, in Figure 
5 we plot the signal-to-noise ratios of (A,E,T) as well as the optimal signal-to-noise ratio. 
For an assumed h = 10 -23 , the SNRs of the three modes are plotted versus frequency. For 
the equal-arm case computed here, the SNRs of A and E are equal across the band. In the 
long wavelength region of the band, modes A and E have SNRs much greater than mode T, 

11 




Iog 10 (frequency, Hz) 



FIG. 3: The LISA Michelson sensitivity curve (SNR = 5) and the sensitivity curve for the optimal 
combination of the data, both as a function of Fourier frequency. The integration time is equal to 
one year, and LISA is assumed to have a nominal armlength L = 16.67sec. 




Iog 10 (frequency, Hz) 

FIG. 4: The optimal SNR divided by the SNR of a single Michelson interferometer, as a function 
of the Fourier frequency /. The sensitivity gain in the low- frequency band is equal to y/2, while it 
can get larger than 2 at selected frequencies in the high-frequency region of the accessible band. 
The integration time has been assumed to be one year, and the proof mass and optical path noise 
spectra are the nominal ones. See the main body of the paper for a quantitative discussion of this 
point. 



where its contribution to the total SNR is negligible. At higher frequencies, however, the T 
combination has SNR greater than or comparable to the other modes and can dominate the 
SNR improvement at selected frequencies. 
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-4 -3.5 -3 -2.5 -2 -1.5 -1 -0.5 

log 10 (frequency, Hz) 

FIG. 5: The SNRs of the three combinations, (A, E, T), and their sum as a function of the Fourier 
frequency /. The SNRs of A and E are equal over the entire frequency band. The SNR of T is 
significantly smaller than the other two in the low part of the frequency band, while is comparable 
to (and at times larger than) the SNR of the other two in the high-frequency region. See text for 
a complete discussion. 

IV. CONCLUSIONS 

The use of Time-Delay Interferometry has shown that LISA has the capability to simulta- 
neously observe gravitational waves in the millihertz band with several, and rather different, 
interferometric data combinations. In this paper we have identified, for a given (but oth- 
erwise arbitrary) gravitational wave signal, the particular interferometric combination that 
gives maximum signal-to-noise ratio. In this context we have actually shown that LISA 
should no longer be regarded as a single-instrument mission, but rather as a network of in- 
terferometer detectors of gravitational radiation working in coincidence. We have identified 
the general expression of the optimal combination of the generators (a,/3, 7), which should 
be used when observing a specified gravitational wave signal. 

Under rather general assumptions on the properties of the noise correlation matrix, and 
for sinusoidal gravitational wave signals that are randomly distributed over the celestial 
sphere and over the polarization states, we have found that the sensitivity gain of the 
optimal data combination over that of the Michelson combination can be significant, and it 
varies over the frequency band accessible by LISA. In the low part of the frequency band 
such improvement is equal to \/2, and it grows to values larger than \f?> at higher frequencies 
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and in small frequency bands centered on frequencies that are integer multiple of the inverse 
of the one-way light-travel time. 

The results derived in this paper will have immediate application to the solution of 
the so called "Inverse Problem" for LISA, that is to say the determination of the source 
location and of the wave's two independent amplitudes from the data LISA will be able to 
generate. We will estimate the accuracies in the determination of the source location and of 
the wave's amplitudes that the optimal combination of the three interferometric responses 
(a, P, 7) derived in this paper will imply. This work is in progress, and will be presented in 
a follow-up publication. 
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